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The  objective  of  this  work  is  to  investigate  the  transport  phenomena  and  performance  of  a  plate  steam 
methanol  micro-reformer.  Micro  channels  of  various  height  and  width  ratios  are  numerically  analyzed  to 
understand  their  effects  on  the  reactant  gas  transport  characteristics  and  micro-reformer  performance. 
In  addition,  influences  of  Reynolds  number  and  geometric  size  of  micro  channel  on  methanol  conversion 
of  micro-reformer  and  gas  transport  phenomena  are  also  explored.  The  predicted  results  demonstrated 
that  better  performance  is  noted  for  a  micro  channel  reformer  with  lower  aspect-ratio  micro  channel. 
This  is  due  to  the  larger  the  chemical  reaction  surface  area  for  a  lower  aspect-ratio  channel  reformer.  It  is 
also  found  that  the  methanol  conversion  decreases  with  increasing  Reynolds  number  Re.  The  results  also 
indicate  that  the  smaller  micro  channel  size  experiences  a  better  methanol  conversion.  This  is  due  to  the 
fact  that  a  smaller  micro  channel  has  a  much  more  uniform  temperature  distribution,  which  in  turn,  fuel 
utilization  efficiency  is  improved  for  a  smaller  micro  channel  reformer. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Micro  fuel  cells  are  widely  regarded  as  the  most  promising 
energy  storage  devices  for  mobiles,  laptops,  and  personal  digital 
assistants  (PDA)  in  the  21st  century  due  to  their  properties  of  high 
energy  density,  low  noise,  and  low  pollution.  Kundu  et  al.  [1]  pre¬ 
sented  that  the  direct  methanol  fuel  cell  (DMFC)  has  still  issues 
due  to  a  low  rate  of  oxidation  and  a  high  crossover  rate.  Therefore, 
methanol  steam  micro-reformers  are  being  developed  for  use  with 
small  proton  exchange  membrane  fuel  cells  (PEMFCs)  to  overcome 
these  problems,  and  could  soon  become  the  new  choice  for  portable 
fuel  cell  applications. 

Flolladay  et  al.  [2]  showed  that  steam  reforming  produces  higher 
yields  of  hydrogen  than  autothermal  reforming  and  partial  oxida¬ 
tion  of  hydrocarbon  fuels,  and  Kolb  et  al.  [3]  demonstrated  that 
the  plate  reactors  have  better  performance  than  cylindrical  reac¬ 
tors  due  to  better  heat  and  mass  transfer.  Therefore,  more  studies 
on  plate  methanol  steam  micro-reformers  are  currently  in  progress. 
As  a  result  of  the  steam  reforming  being  an  endothermic  reaction, 
several  researchers  have  used  electrical  heat  sources  to  supply  heat 
flux  to  steam  micro-reformers  [4-6].  Ha  et  al.  [4]  fabricated  a  PDMS 
(poly  dimethylsiloxane)  micro-reformer,  and  used  a  heater  to  sup¬ 
ply  heat  flux.  Their  experimental  results  showed  that  methanol 
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conversion  was  about  30-40%  in  the  operating  temperatures  of 
180-240  °C.  The  methanol  conversion  and  hydrogen  production  of 
a  plate-type  micro-reformer  with  micro  channels  patterned  on  the 
micro-reformer  was  examined  experimentally  by  Lim  et  al.  [5].  The 
results  showed  approximately  78%  conversion  of  methanol,  with  a 
hydrogen  production  rate  of  3L/h.  Kwon  et  al.  [6]  used  silicon  to 
fabricate  a  reformer  and  preferential  oxidizer  (PrOx)  reactor.  The 
reformer  produced  hydrogen  to  supply  the  fuel  cell.  The  measured 
data  disclosed  that  the  reformer  generated  27  cm3 /min  of  hydro¬ 
gen  and  that  the  CO  was  totally  removed  from  the  gas  by  the  PrOx 
device.  When  the  fuel  cell  was  operating  at  0.6  V,  the  power  density 
was  230  mW/cm2. 

More  recently,  several  studies  have  used  catalytic  burners  to 
supply  thermal  energy  to  the  entire  micro-reformer  [7-9].  Kwon 
et  al.  [7]  utilized  silicon  fabrication  technology  to  set  up  the 
reformer  and  the  catalytic  burner.  Their  results  indicated  that 
the  methanol  combustion  reaction  with  the  catalytic  burner  suc¬ 
cessfully  generated  heat  to  maintain  reformer  temperature.  The 
methanol  steam  reformer  combined  with  the  catalytic  burner  pro¬ 
duced  73%  hydrogen  with  65%  conversion  of  the  methanol.  A  micro 
channel  methanol  steam  reformer  which  included  a  vaporizer,  heat 
exchanger,  catalytic  combustor  and  steam  reformer  to  produce 
hydrogen  for  a  PEMFC  was  conducted  experimentally  by  Park  et 
al.  [8].  Their  results  showed  that  at  a  temperature  of  250  °C,  the 
reformer  could  produce  a  gas  flow  rate  of  450  ml/min,  with  a  gas 
composition  of  73.3%  H2  at  the  micro  channel  outlet.  Yoshida  et  al. 
[9]  fabricated  a  micro  channel  methanol  reformer  integrated  with  a 
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Nomenclature 

Q  concentration  of  species  i  (mol  m-3 ) 

Dk  mass  diffusion  coefficient  (m2  s-1 ) 

Dp  catalyst  particle  diameter  (m) 

Ea  activation  energy  (J  mol-1 ) 

H  micro  channel  height  (m) 

hi  species  enthalpy 

/J,  I<  grid  points  in  the  x,  y  and  z  directions,  respectively 
kp  permeability  (m2) 

/<i  pre-exponential  factor  for  steam  reforming 

k 2  pre-exponential  factor  for  the  reverse  water  gas  shift 

L  micro  channel  length  (m) 

Mf  mole  fraction  of  species  i 

MWI  molecular  weight  of  species  i  (g  mol-1 ) 
rrii  mass  fraction  of  species  i 

Nr  number  of  chemical  species  in  the  reaction 

p  pressure  (Pa) 

qHl  hydrogen  production  rate  at  outlet  (cm3  min-1 ) 

R  universal  gas  constant 

T  temperature  (°C) 

Tw  wall  temperature  (°C) 

W  micro  channel  width  (m) 

x,  y,  z  coordinates  (m) 

u,  v ,  w  velocity  components  in  the  x,  y  and  z  directions, 

respectively,  (ms-1) 

Greek  symbols 

a  thermal  diffusivity  (m2  s-1 ) 

y  aspect  ratio  (height  and  width  ratio,  HjW) 

<5i  catalyst  layer  height  (m) 

82  flow  channel  height  (m) 

s  porosity 

p  methanol  conversion 

fi  viscosity  (Ns m-2) 

p  density  (kg  m-3) 

Subscripts 

u  x-direction 

v  y-direction 

w  z-direction 

0  inlet 


combustor  and  a  micro  channel  vaporizer.  An  appropriate  depth  of 
the  micro  channel  vaporizer  could  generate  high  yields  of  hydrogen 
was  found  in  their  work. 

There  are  many  research  works  about  the  developed  catalysts 
used  in  the  reforming  reaction  using  different  kinetics  of  the 
methanol  reforming  reactions.  Peppley  et  al.  [10]  studied  the  reac¬ 
tion  network  for  methanol  steam  reforming  over  a  Cu/Zn0/Al203 
catalyst  form  BASF  K3-110.  The  steam  reforming  of  methanol  over 
a  Cu/Zr02/Ce03  catalyst  was  investigated  by  Mastalir  et  al.  [11]. 
The  kinetic  model  suggested  for  the  transformation  involved  the 
reverse  water-gas  shift  and  methanol  decomposition,  in  addition 
to  the  steam  reforming  methanol  reaction.  Lee  et  al.  [12]  carried 
out  a  kinetic  study  of  methanol  steam  reforming  over  a  commercial 
catalyst  Cu0/Zn0/Al203. 

In  order  to  reduce  the  research  cost  and  design  duration, 
modeling  and  simulation  are  often  being  used  to  obtain  a  bet¬ 
ter  understanding  of  the  processes  in  a  methanol  reformer  at 
present.  Known  from  the  literature  survey,  several  numerical  mod¬ 
els  have  been  presented,  with  the  simplest  methods  using  a 
one-dimensional  model  to  describe  the  methanol  reformer  con¬ 
version  and  the  heat  and  mass  transport  phenomena  [13-18]. 


Kawamura  et  al.  [13]  proposed  a  micro  channel  model  using  mass 
and  heat  balances  to  analyze  the  transport  phenomena  in  a  plate 
reformer.  They  successfully  simulated  methanol  conversion  and  gas 
concentration  distributions  along  the  micro  channel.  A  mathemat¬ 
ical  model  of  a  plate  methanol  steam  reformer  to  investigate  the 
transport  phenomena  was  studied  numerically  by  Kim  and  Kwon 
[14].  Their  results  showed  that  smaller  reformer  volumes  required 
a  higher  heat  flux.  Varesano  et  al.  [15]  used  a  one-dimensional 
transient  mathematical  model  to  study  the  transport  behaviors  in 
a  steam  reforming  reformer  with  a  burner.  The  transient  charac¬ 
teristics  of  the  reformer  were  examined  by  Amphlett  et  al.  [19] 
using  the  reforming  kinetics.  Pattekar  and  Kothare  [16]  developed 
one-dimensional  models  for  radial  and  micro  channel  methanol 
reformers.  Their  results  demonstrated  that  radial  reformers  have 
better  hydrogen  production  rates  and  lower  pressure  drops  than 
micro  channel  reformers.  A  transient  model  for  packed-bed  cat¬ 
alytic  methanol  reformers  was  performed  by  Stamps  and  Gatzke 
[17].  They  found  that  agreement  of  the  theoretical  and  experimen¬ 
tal  results  for  the  temperature,  flow  rate  and  CO  concentration  data. 
Lattner  and  Harold  [  18  ]  have  numerically  and  experimentally  inves¬ 
tigated  a  bench-scale  fixed-bed  methanol  reformer  for  autothermal 
reforming.  This  system  was  also  used  by  Pepply  et  al.  [20]  and  Reitz 
et  al.  [21]  who  presented  kinetic  rate  expressions  for  the  simula¬ 
tions.  Their  results  showed  agreement  between  the  theoretical  and 
experimental  results. 

As  for  the  two-dimensional  simulation  about  the  methanol 
reformers,  the  literature  survey  has  existed  in  the  past  decade.  Suh 
et  al.  [22]  proposed  a  cylindrical  mathematical  model  of  a  packed 
bed  reformer  to  investigate  the  heat  and  mass  transfer  for  differ¬ 
ent  species  in  a  reformer.  The  results  indicated  that  the  predicted 
methanol  conversion  rates  were  compared  well  with  the  exper¬ 
imental  data.  A  mathematical  model  of  a  cylindrical  reformer  to 
investigate  the  methanol  conversion  in  a  packed  bed  reformer  and 
a  wall-coated  reformer  was  proposed  by  Karim  et  al.  [23,24].  Their 
results  showed  that  the  minimum  reformer  diameter  yielded  the 
highest  catalyst  activity  and  smallest  temperature  gradient.  The 
results  also  showed  that  wall-coated  reformers  have  higher  cat¬ 
alytic  activity  and  lower  pressure  drops  than  packed-bed  reformers 
do. 

There  are  many  three-dimensional  simulations  of  methanol 
reformers  in  the  literature.  A  cylindrical  model  of  the  reformer 
which  comprised  of  a  methanol  steam  reformer  and  a  CO  methana- 
tor  to  simulate  the  conversion  and  temperature  distributions  in 
the  reformer  was  investigated  by  Cao  et  al.  [25].  Pan  and  Wang 
[26]  developed  a  numerical  model  for  a  plate-fin  reformer  which 
integrated  endothermic  and  exothermic  reactions  into  one  unit. 
Their  numerical  model  accurately  predicted  the  methanol  conver¬ 
sion  and  the  gas  distributions.  Numerical  simulation  of  transport 
phenomena  of  the  plate  micro-scale  methanol  reformer  was  exam¬ 
ined  by  Pattekar  and  Kothare  [27].  Park  et  al.  [28]  developed 
three-dimensional,  quasi  three-dimensional  and  one-dimensional 
models  to  study  the  plate  methanol  reformers.  Cao  et  al.  [29] 
presented  kinetic  rate  expressions  and  developed  a  homogenous 
model  of  a  micro  channel  reformer  to  simulate  the  temperature 
distributions  in  a  micro  channel. 

The  geometric  design  of  a  reactor  is  one  of  the  most  important 
issues.  Appropriate  reactor  geometry  can  improve  the  reactant  gas 
transport  and  the  efficiency  of  thermal  management  [20,21,30-36]. 
To  this  end,  Kundu  et  al.  [37]  used  different  flow  configurations 
including  serpentine  and  parallel  flow  fields  to  improve  plate 
methanol  reformer  performance.  As  stated  above,  while  many  stud¬ 
ies  have  investigated  the  effects  of  reactor  radius  on  cylindrical 
reactor  performance  [23,24,32-36],  few  studies  have  reported  on 
the  flow  channel  designs  of  plate  methanol  steam  micro-reformers. 
Therefore,  based  on  flow  channel  designs,  various  aspect  ratios 
of  micro  channels  on  plate  methanol  steam  micro-reformers  can 
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potentially  enhance  fuel  utilization.  In  this  work,  flow  channels 
with  various  aspect  ratios  (height  and  width  ratios)  and  geomet¬ 
ric  size  are  numerically  examined  the  transport  phenomena  in  a 
micro  channel  reformer.  In  addition,  the  thermo-fluid  parameters 
(Reynolds  number  and  wall  temperature)  are  also  investigated  to 
examine  their  effects  on  the  methanol  conversion  and  efficiency  of 
micro  channel  reformers. 

2.  Analysis 

In  the  present  study,  a  three-dimensional  micro  channel  model 
of  the  plate  methanol  steam  micro-reformer  was  proposed  to 
analyze  the  local  transport  phenomena  and  micro-reformer  per¬ 
formance.  The  micro  channel  is  comprised  of  the  flow  channel  and 
a  catalyst  layer.  Therefore,  the  governing  equations  include  mass, 
momentum,  energy  and  species  equations.  For  simplicity,  the  fol¬ 
lowing  assumptions  are  made: 

(1 )  The  flow  is  steady  state. 

(2)  The  inlet  fuel  is  an  ideal  gas. 

(3)  The  flow  is  laminar  and  incompressible. 

(4)  The  catalyst  layer  is  isotropic. 

(5)  The  chemical  reaction  occurs  only  in  the  catalyst  layer. 

(6)  Thermal  radiation  and  conduction  in  the  gas  phase  are  negligi¬ 
ble  compared  to  convection. 

To  reduce  the  computing  time,  the  symmetric  micro  channel  is 
considered  only  in  this  work.  The  schematic  diagram  of  this  work 
is  shown  in  Fig.  1.  The  various  cases  for  different  aspect  ratios  are 
shown  in  Table  1.  The  aspect  ratios,  y ,  are  defined  as  follows: 


Table  1 

The  cases  with  various  aspect-ratio  micro  channels  used  in  this  work. 


where  H  and  W  are  the  micro  channel  height  and  width,  respec¬ 
tively.  The  corresponding  hydraulic  diameters  in  Table  1  are  fixed 
to  be  0.286  mm. 

The  governing  equations  under  the  assumptions  stated  above 
can  be  defined  as  follows: 

Continuity  equation: 


ep 


( du  dv  9w\ 
y  dx  +  dy  +  dz  J 


=  0 


(2) 


X-momentum  equation: 


du  du  du 


JP 


£P  ~  ~  £*r+£ll  (3) 


dx  dy  dz 

Y-momentum  equation: 


'  dx 


d2u  d2u  d2u 


f  dv  dv  dv\  dp  .  „  v  „  v  „  „  x 

sp  [ute+vfy+wte ) -~%+£lx  +  +  1  +s"  W 

Z-momentum  equation: 


dx2  dy 2  dz2 


d2v  d2v  d2v 


dw  dw 


dw 


splulK+vW  +  wte 


dp  f  d2W  d2W  d2w\ 

=  _e_+e^_  +  _  +  _j 


+  5w 


(5) 


SU,SV  and  Sw  denote  source  terms  based  on  the  Darcy  drag  forces 
in  the  x,  y  and  z  directions  imposed  by  the  pore  walls  on  the  fluid 
which  usually  cause  significant  pressure  drops  across  the  porous 
media.  In  this  work,  the  flow  channel  is  not  a  porous  medium,  there¬ 
fore,  S U,SV  and  Sw  are  zero  in  the  flow  channel.  Flowever,  the  catalyst 
layer  is  a  porous  medium,  so  the  source  terms  in  the  momentum 
equations  for  the  fluid  flow  in  the  porous  media  are: 


lip  2 

(6) 

Sv  =  -  J^-  -  a/u2  +  1/2  +  w2 

lip  2 

(7) 

s„.-'^-^y„w+w2 

lip  2 

(8) 

where  kp  is  the  permeability  and  j3  is  inertial  loss  coefficient  in  each 

component  direction  [38]. 

,  Dp2e3 

kp  — 

0) 

150(1 -e)2 

«  3.5(1  -e) 

DpS3 

(10) 

Species  equation: 


dQ  dQ  dq 

U^-L  +w^ 

dx  dy  dz 


=  De\'f 


fd2q  d2Cj  d2ct  \ 

Y  dx2  dy2  dz2  J 


+  sSc 


(11) 


In  the  species  equation,  Q  denotes  the  concentration  of  the  ith 
species:  the  calculations  included  CFI3OFI,  FI20,  H2,  C02  and  CO.  In 
Eq.  ( 11 ),  the  effective  mass  diffusivity  Deff  is  expressed  as  Deff  =  sTDk, 
where  £  is  the  porosity  of  the  medium  and  r  is  the  tortuosity  of  the 
porous  medium.  In  this  work,  r  is  set  to  be  1,  and  the  porosity  s  is 
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expressed  as  0.38  and  1.00,  respectively,  in  the  catalyst  layer  and  the 
flow  channel.  In  the  present  study,  there  is  no  chemical  reaction  in 
the  flow  channel,  therefore,  Sc  is  zero  in  flow  channel.  In  the  catalyst 
layer,  the  source  term  of  the  concentration  equation,  Sc,  stands  for 
the  modified  concentration  term  produced  by  the  chemical  reac¬ 
tion.  Sc  is  the  sum  of  the  Arrhenius  reaction  sources  over  the  NR 
reactions  that  the  species  participle  in: 

Nr 

Sc=MwJ  J2  RUr  (12) 

r=SR,rWSG 

where  Mwl  is  the  molecular  weight  of  species  i  and  Ri<r  is  the 
Arrhenius  molar  rate  of  creation  and  destruction  of  species  i  in  the 
reaction  r.  According  to  the  simple  reaction  model  of  Purnama  et  al. 
[39],  the  steam  reforming  reaction  is  much  faster  than  the  decom¬ 
position  and  water-gas  shift  reaction.  Therefore,  only  the  steam 
reforming  reaction,  Eq.  (13)  and  the  water-gas  shift  reaction,  Eq. 
(14)  are  considered  in  this  work. 

CH3OH  +  H20-^C02  +  3H2  (13) 

C0  +  H20^C02+H2  (14) 

k-2 

Using  the  model  for  steam  reforming  of  methanol  by  Pan 
and  Wang  [26],  who  presented  the  following  mole  rate  of  cre¬ 
ation/destruction  of  species  i  in  the  SR  reaction  and  the  rWGS 
reaction: 

RU  SR  =  ^1  Cm*  OH  ^H240  eXP  ( ' “  ^  5 ) 

Ri, rWSG  =  ^2Cco2Ch2  exp  “  ^-2CcqCh20  exp  O6) 


Table  2 

Parameters  used  in  this  study. 


Micro  channel  length  L  (m)  [40] 

Micro  channel  width  W (m)  [40] 

Micro  channel  height  H  (m)  [40] 

Catalyst  layer  thickness  <5i  (m) 

Flow  channel  height  82  (m) 

Average  inlet  temperature  (°C)  [40] 

Operating  pressure  (atm)  [40] 

Activation  energy  for  steam  reforming  (J  mol-1 )  [26] 

Activation  energy  for  reverse  water  gas  shift  (J  mol-1 )  [26] 
Catalyst  density  (kgnrr3)  [23] 

Catalyst  layer  porosity  [27] 

Catalyst  permeability  (m2)  [27] 

3.3  x  10-2 

5.0  x  10-4 

2.0  x  10-4 

3.0  x  10-5 

1.7  x  10-4 

120 

1 

7.6  x  104 

1.08  x  105 

890 

0.38 

2.379  x  10-12 

C? 

II 

cr 

(21) 

11 

(22) 

At  the  micro  channel  outlet  (x  =  L),  fully  developed  flow  is 
assumed  and  the  boundary  conditions  of  velocity,  temperature  and 
concentration  can  be  expressed  as: 

du  dC  dT  „ 

dx  dx  dx 

(23) 

At  the  interface  between  the  flow  channel  and  the  insulated 
plates  (y  =  0),  the  velocities  and  the  temperature  and  concentration 
gradients  are  zero. 

dCi  dT 

U=V=W  =  -^4.  =  _  =  0 

dy  dy 

(24) 

At  the  interface  between  the  flow  channel  and  the  catalyst  layer 
(y  =  Hi  -  82 ),  the  velocities,  temperature,  species  concentration,  and 
species  flux  are  continuous. 


Where  the  steam  reforming  reaction  is  a  non-reversible  reaction 
and  the  reverse  water-gas  shift  reaction  is  reversible.  The  con¬ 
stants,  /<i  and  k2,  are  forward  rate  constants  for  the  steam  reforming 
reaction  and  the  reverse  water-gas  shift  reaction,  respectively.  The 
constant,  /<_ 2,  is  the  backward  rate  constant  for  the  water-gas  shift 
reaction. 

Energy  equation: 


.  3T  dT  37  \  ,  f  d2T  d2T  32I  .  „ 

pCp  l  u  &  +  vty  + w  &  )  -  keff  (  W  +  w  +  a?  )  +St 


(17) 
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Ci,y=(Hj-S2)+  ~  Ci,y=(Hj-S2r 


D. 


eff ,(H-S2f 


9Q 
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In  the  energy  equation,  the  effective  thermal  conductivity  in  the 
porous  medium,  keff,  is  the  volume  average  of  the  fluid  conductivity 
and  solid  conductivity. 


“  ^eff,(H-52r 


3Q 

3y 


y=(H-S2r 


(27) 


keff  =  skf  +  (1  —  s)1<s  (18) 

where  kf  is  fluid  phase  thermal  conductivity,  ks  is  the  solid  medium 
thermal  conductivity  and  s  is  the  porosity  of  the  medium. 

Source  term,  St,  in  Eq.  (17)  includes  the  heat  flux  due  to  the 
chemical  reactions.  In  this  work,  the  catalyst  layer  experiences 
exothermic  chemical  reactions,  therefore,  St  for  the  chemical  reac¬ 
tions  is: 


(19) 


T  _  T  /  dT 

ly=(H^-h)+  ~  y=(H-[  —82)~  ’  Keff,(H-82)+  fy 


y=(H-82 )+ 


=  k 


8T 

eff ,y=(H-S2r  ^ y 


y=(H-S2r 


(28) 


At  the  interface  between  the  heated  wall  and  the  catalyst  layers, 
the  velocities  and  the  concentration  gradient  are  assumed  to  be 
zero,  and  the  temperature  is  assumed  to  be  equal  to  the  catalyst 
temperature: 


where  is  the  enthalpy  of  formation  of  species  i,  and  R,  is  the 
volumetric  rate  of  creation  of  species  i. 

The  boundary  conditions  of  the  present  computation  include 
those  at  the  inlet,  outlet,  wall,  and  the  interface  between  the  flow 
channel  and  the  catalyst  layer.  The  inlet  flow  velocity  is  constant, 
the  inlet  gas  compositions  are  constant,  and  the  inlet  temperature 
is  constant. 

(20) 


U  =  v  =  w=—L=0,  T  =  TW  (29) 

3y 

In  the  present  study,  the  mass  fractions  of  the  inlet  reactant  gas 
including  methanol  vapor  and  water  vapor  of  0.38  and  0.62,  respec¬ 
tively,  were  tested.  Thus,  the  molar  ratio  of  H20/CH30H  was  kept 
constant  at  1.1.  The  inlet  flow  velocity  of  0.266  m/s  used  and  there¬ 
fore,  the  corresponding  Reynolds  number  Re  was  10.98  for  each  test. 
The  physical  properties  of  the  micro  channels  are  listed  in  Table  2. 
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Table  3 

Mole  fractions  of  methanol  for  the  various  grid  tests  at  different  axial  locations. 


IxJxK  x(m) 


0.005 

0.010 

0.015 

0.020 

0.025 

0.030 

41  x  18  x  18 

0.418 

0.372 

0.334 

0.303 

0.276 

0.253 

51  x  25  x  23 

0.419 

0.372 

0.335 

0.304 

0.277 

0.254 

71  x  35  x  33 

0.419 

0.373 

0.335 

0.303 

0.277 

0.253 

3.  Numerical  method 


A  generalized  form  of  the  transport  equation  for  mass,  momen¬ 
tum,  energy  can  be  expressed  in  a  conservative  form  as  follows: 

^ JpcpdV  =  Jp<pVdA  =  j  rv<pdA  +  Js^dV  (30) 


where  t  is  the  time,  0  is  a  general  dependent  variable,  5^  is 
the  source  per  unit  volume  and  p  is  the  density.  All  equations 
were  numerically  solved  using  the  commercial  fluid  dynamics  pro¬ 
gram,  Fluent™.  The  SIMPLE  algorithm  was  employed  to  solve  the 
convection-diffusion  equations. 

The  grid  independence  was  examined  in  preliminary  test  runs. 
Three  grid  configurations  were  evaluated  for  the  micro  channel  of 
the  plate  methanol  steam  micro-reformer  at  a  wall  temperature  of 
200  °C.  The  numbers  of  grid  lines  in  the  x,  y  and  z  directions  were: 
41  x  18  x  18,  51  x  25  x  23,  and  71  x  35  x  33.  The  influence  of  grid 
lines  on  the  local  methanol  mole  fraction  is  shown  in  Table  3.  The 
deviations  of  methanol  mole  friction  are  0.4%  for  grids  41  x  18  x  18 
and  51  x  25  x  23,  and  0.4%  for  grids  51  x  25  x  23  and  71  x  35  x  33. 
Grid  51  x  25  x  23  was,  therefore,  chosen  for  the  simulation  in  the 
present  study  as  a  tradeoff  between  accuracy  and  CPU  computation 
time.  The  convergence  criterions  for  the  normalized  residuals  for 
each  variable  were  restricted  to  less  than  10-6. 

The  validation  of  the  numerical  results  is  performed  by  com¬ 
paring  the  present  predictions  with  previous  experimental  results. 
The  geometric  dimensions  and  physical  properties  of  the  micro 
channel  are  listed  in  Table  2.  Fig.  2  shows  the  comparison  of  the 
present  predictions  and  experimental  results.  The  solid  symbols 
denote  the  experimental  results  of  Park  et  al.  [40]  and  the  curve  is 
the  present  prediction.  The  results  show  that  the  numerical  results 
agree  reasonably  well  with  the  experimental  data. 

4.  Results  and  discussion 

In  this  work,  the  numerical  results  are  obtained  for  a  micro  chan¬ 
nel  of  a  plate  methanol  steam  micro-reformer.  Using  previously 
stated  the  numerical  model,  the  effects  of  aspect  ratios  of  micro 


Fig.  3.  Variations  of  the  mole  fractions  of  the  various  species  along  the  micro  channel 
center  line  (y  =  0.5). 


channel  on  methanol  conversion  and  transport  phenomena  were 
emphasized.  For  aspect  ratio  of  y  =  0.5,  Fig.  3  presents  the  local  dis¬ 
tributions  of  the  different  species  at  wall  temperatures  of  200  and 
260  °C  along  the  center  of  the  micro  channel.  Overall  inspection  of 
Fig.  3  disclosed  that  the  CH3OH  and  H20  mole  fractions  decrease 
along  the  micro  channel,  while  the  H2,  C02  and  CO  mole  fractions 
increase  along  the  center  of  the  micro  channel.  The  results  demon¬ 
strate  that  as  the  mole  fractions  of  the  products  increase  as  the  wall 
temperature  increases.  The  results  also  show  that  the  methanol 
conversion  is  about  49%  for  a  wall  temperature  of  200 °C,  with  a 
gas  composition  of  24%  CH3OH,  28%  H20,  36%  H2,  and  12%  C02  at 
the  micro  channel  outlet.  However,  the  CO  concentration  is  only 
244  ppm.  For  a  wall  temperature  of  260  °C,  the  results  indicates 
that  the  methanol  conversion  is  greater  than  96%,  with  a  product 
gas  composition  of 74.05%  H2 , 24.28%  C02 ,  and  1.67%  CO  at  the  micro 
channel  outlet.  For  the  PEMFC,  the  CO  concentration  must  be  less 
than  10  ppm  which  can  be  achieved  using  a  CO  oxidation  reactor. 
The  utilization  of  the  preferential  oxidizer  (PrOx)  reactor  or  water- 
gas-shift  reaction  can  reduce  the  CO  concentration  in  the  gas  from 
the  methanol  micro-reformer.  Kwon  et  al.  [6]  used  a  reformer  and 
PrOx  reactor,  and  showed  that  the  CO  produced  was  totally  removed 
from  the  gas  by  the  PrOx  device. 

Fig.  4  presents  the  effects  of  the  aspect  ratios  of  micro  channel 
on  methanol  conversion  and  CO  concentration  (ppm).  In  this  work, 
the  catalyst  thicknesses  are  fixed  to  be  30  |jim.  A  careful  examina¬ 
tion  of  Fig.  4  reveals  that  the  methanol  conversion  increases  with 
an  increase  in  wall  temperature.  Thus,  the  methanol  conversion 
can  be  improved  by  increasing  wall  temperature,  which  in  turn, 
increases  the  chemical  reaction  rate.  It  is  also  found  form  Fig.  4  that 
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Fig.  2.  Comparison  of  predicted  methanol  conversion  with  previous  experimental 
data  of  Park  et  al.  [40]. 


Fig.  4.  Effects  of  aspect  ratios  of  micro  channel  and  wall  temperature  on  methanol 
conversion  and  CO  concentration  (ppm)  at  outlet  of  micro  channel. 
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Fig.  5.  Effects  of  aspect  ratios  of  micro  channel  on  the  local  methanol  conversion 
along  center  line  of  micro  channel  at  Tw  =  200  °C  and  Tw  =  260  °C. 


the  methanol  conversion  increases  with  a  decrease  in  aspect  ratio 
due  to  the  large  chemical  reaction  area  for  a  low  aspect-ratio  micro 
channel.  This  implies  that  better  performance  is  noted  for  a  lower 
aspect-ratio  micro  channel  reformer.  The  CO  concentration  (ppm) 
leaving  the  micro  channel  was  also  studied.  In  Fig.  4,  the  results 
show  that  the  CO  concentration  increases  with  increasing  wall  tem¬ 
perature.  This  can  be  made  plausible  by  noting  the  fact  that  the 
reaction  rate  of  the  endothermic  reverse  water-gas-shift  reaction 
increases  as  the  wall  temperature  increases.  Additionally,  it  is  also 
found  that  the  CO  concentration  (ppm)  increases  with  decreasing 
aspect  ratios.  A  height/width  ratio  of  y  =  0.25  and  a  wall  tempera¬ 
ture  of  260  °C  yielded  the  better  methanol  conversion  (greater  than 
98%),  but  the  CO  concentration  was  also  higher,  and  was  greater 
than  16,000  ppm.  Hence,  the  CO  concentration  in  the  outlet  gases 
must  be  reduced  for  further  use  in  a  PEMFC. 

The  effects  of  aspect  ratios  on  local  methanol  conversion  dis¬ 
tributions  along  the  micro  channel  center  line  at  TW  =  200°C  and 
Tw  =  260°C  are  presented  in  Fig.  5.  An  overall  inspection  of  Fig.  5 
reveals  that  the  aspect  ratios  of  micro  channel  have  a  consider¬ 
able  impact  on  the  local  methanol  conversion  distributions  along 
the  center  line.  It  is  found  that  the  predicted  methanol  conversion 
increases  along  the  micro  channel.  In  addition,  methanol  conver¬ 
sion  efficiency  increases  with  decreasing  aspect  ratios  of  micro 
channel.  As  for  the  effects  of  wall  temperature,  better  methanol 
conversion  is  noted  for  a  case  with  a  higher  wall  temperature  owing 
to  a  stronger  chemical  reaction  rate. 

Fig.  6  shows  the  effects  of  the  aspect  ratios  of  micro  channel 
on  the  distribution  of  the  H2  and  CO  mole  fractions,  respectively, 
along  the  micro  channel  center  line.  It  is  clear  in  Fig.  6(a)  that  the 
H2  mole  fraction  along  the  micro  channel  center  line  represents  the 
methanol  conversion,  with  a  higher  H2  fraction  indicating  higher 
methanol  conversion  rates.  Thus,  the  variation  of  the  H2  fraction  is 
opposite  to  that  of  the  CH3OH  mole  fraction.  It  is  clearly  observed 
from  Fig.  6(a)  the  H2  mole  fraction  increases  as  the  aspect  ratio 
decreases  and  wall  temperature  increases.  The  trends  of  the  vari¬ 
ations  in  Fig.  6(b)  can  be  interpreted  in  a  similar  way  to  the  data 
in  Fig.  5.  This  is  because  that  a  higher  methanol  conversion  results 
in  a  higher  CO  production.  These  phenomena  are  more  obvious  at 
a  higher  wall  temperature  (260  °C),  because  the  methanol  micro¬ 
reformer  produces  more  CO  for  a  higher  wall  temperature. 

Fig.  7(a)  presents  the  reactant  gas  velocity  distributions  in  the 
center  of  the  micro  channel  along  the  direction  of  flow  from  the 
inlet  to  the  outlet  for  wall  temperatures  of  200  and  260  °C  for  the 
various  aspect  ratios.  It  is  clearly  seen  in  Fig.  7(a)  that  the  velocity 
slowly  increases  from  the  inlet  to  the  outlet.  The  wall  tempera¬ 
ture  of  260  °C  generates  higher  gas  velocity  distributions  than  that 
of  200  °C.  The  results  show  that  the  methanol  conversion  of  the 


micro-reformer  is  independent  of  the  velocity  distributions  in  the 
center  of  the  micro  channel  due  to  the  chemical  reaction  dominat¬ 
ing  methanol  conversion.  Then,  to  explore  the  pressure  loss  caused 
by  the  various  aspect  ratios,  the  local  pressure  losses  (the  differ¬ 
ence  between  local  and  inlet  pressures)  along  the  micro  channel 
center  line  for  wall  temperatures  of  200  and  260  °C  under  various 
aspect  ratios  are  presented  in  Fig.  7(b).  In  Fig.  7(b),  it  is  found  that 
the  local  pressure  loss  increases  along  the  flow  channel.  A  careful 
examination  of  Fig.  7(b)  discloses  that  larger  pressure  loss  is  noticed 
for  a  case  with  a  higher  wall  temperature  (260  °C).  Additionally,  the 
results  show  that  larger  pressure  loss  is  found  for  a  smaller  aspect 
ratio  (y  =  0.25).  The  higher  pressure  losses  mean  that  excess  work 
has  to  be  done  in  pushing  the  gases  through.  Therefore,  reducing 
pressure  losses  is  one  of  the  most  important  issues. 

Effects  of  the  aspect  ratios  on  the  temperature  distributions  are 
shown  in  Fig.  8(a).  It  is  clearly  seen  that  the  fluid  temperature 
increases  along  the  channel  as  a  consequence  of  chemical  reaction. 
When  x/L  >  0.2,  the  fluid  temperature  increase  insignificantly  along 
channel  because  there  is  a  thermal  equilibrium  beyond  this  point. 
It  is  also  shown  that  the  temperature  rises  with  the  decrease  in 
aspect  ratios,  y.  As  the  y  decreases,  the  temperature  distributions 
become  much  more  uniform.  Therefore,  a  smaller  aspect  ratio  leads 
to  a  better  methanol  conversion.  Afterwards,  to  explore  the  change 
in  steam  reforming  reaction  rate  caused  by  the  various  aspect  ratios 
of  the  micro  channels,  the  local  steam  reforming  reaction  rates 
along  the  interface  between  the  catalyst  layer  and  the  flow  chan¬ 
nel  are  presented  in  Fig.  8(b).  For  a  wall  temperature  of  260  °C,  the 
steam  reforming  reaction  rate  first  rises  rapidly,  and  then  declines 
steadily  after  x/L  =  0.001.  For  x/L  <  0.05,  the  steam  reforming  reaction 
rate  increases  with  decreasing  aspect  ratios.  For  x/L  >0.05,  larger 
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Fig.  6.  Effects  of  aspect  ratios  of  micro  channel  at  7W  =  200  °C  and  TW  =  260°C  on 
(a)  local  H2  mole  fraction  and  (b)  local  CO  mole  fraction  along  center  line  of  micro 
channel. 
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Fig.  7.  Effects  of  aspect  ratios  of  micro  channel  at  Tw  =  200  °C  and  Tw  =  260  °C  on  (a) 
local  velocity  and  (b)  local  pressure  along  center  line  of  micro  channel. 


aspect  ratios  can  improve  the  steam  reforming  reaction  rate,  as  a 
result  of  the  methanol  being  consumed  in  the  forward  region  of 
micro  channel  inlet.  For  a  lower  wall  temperature  of  200  °C,  the 
change  in  steam  reforming  reaction  rate  follows  a  similar  trend  to 
that  of  the  wall  temperature  260  °C.  In  addition,  it  is  found  that 
a  higher  steam  reforming  reaction  rate  is  found  for  a  case  with  a 
higher  wall  temperature  due  to  a  stronger  steam  reforming  reac¬ 
tion. 

Fig.  9  demonstrates  the  effects  of  the  Reynolds  number  on  the 
methanol  conversion  and  the  H2  production  rate  of  the  micro 
channel  in  the  plate  methanol  steam  micro-reformer  for  wall  tem¬ 
peratures  of  260  °C.  An  overall  inspection  of  Fig.  9  discloses  that 
the  methanol  conversion  increases  with  decreasing  Reynolds  num¬ 
ber  Re.  This  can  be  made  plausible  by  noting  the  fact  that  the  fuel 
can  spend  more  time  in  the  micro  channel  for  a  lower  Reynolds 
number  Re.  As  for  the  H2  production  rate,  the  results  show  that 
the  H2  production  rate  increases  as  the  Re  increases.  Additionally, 
the  effects  of  the  aspect  ratios  on  the  methanol  conversion  and  FI2 
production  rate  are  shown  in  Fig.  9.  It  is  clear  that  the  methanol 
conversion  increases  with  a  decrease  in  aspect  ratio.  When  y  =  0.25 
with  a  Re  of  22,  the  methanol  conversion  was  greater  than  86%,  and 
the  H2  production  rate  was  231  cm3 /min.  The  plate  methanol  steam 
micro-reformer  produced  hydrogen  to  supply  fuel  cell,  Park  et  al. 
[40]  presented  that  the  H2  production  rate  was  only  186cm3/min, 
which  could  supply  a  15  W  fuel  cell. 

Effects  of  geometric  size  on  the  transport  phenomena  and  per¬ 
formance  of  micro-reformer  are  important.  To  this  end,  the  effects 
of  micro  channel  size  on  temperature  distributions  are  shown  in 
Fig.  10.  In  Fig.  10,  two  micro  channel  sizes  were  tested:  micro  chan¬ 
nel  (I)  33  mm(L)  x  1.5  mm( IV)  x  0.75  mm(H)  and  micro  channel  (II) 


Fig.  8.  Effects  of  aspect  ratios  of  micro  channel  at  Tw  —  200  °C  and  Tw  —  260  °C  on  (a) 
local  temperature  along  center  line  of  micro  channel  and  (b)  local  reaction  rates  of 
methanol  steam  reforming  along  interface  between  flow  channel  and  catalyst  layer. 

33  mm(L)  x  0.429  mm(W)  x  0.214  mm(H),  respectively.  The  aspect 
ratios  were  all  fixed  to  be  0.5,  with  the  catalyst  thickness  being 
30  |jim.  Different  micro  channel  sizes  have  different  hydraulic  diam¬ 
eters.  Micro  channel  (I)  has  a  greater  hydraulic  diameters  than  micro 
channel  (II).  It  is  found  from  Fig.  10  that  the  temperature  distri¬ 
bution  increases  along  the  micro  channel  as  a  consequence  of  the 
chemical  reaction.  For  a  smaller  micro  channel  size,  the  tempera¬ 
ture  distribution  is  much  more  uniform  due  to  the  shorter  thermal 
entrance  length.  This  kind  of  uniform  temperature  distribution 
improves  the  chemical  reaction  rate.  As  a  result,  the  smaller  micro 
channel  leads  to  better  methanol  conversion.  Comparing  the  local 
temperature  distributions  for  wall  temperatures  of  200  and  260  °C, 
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Fig.  9.  Effects  of  aspect  ratios  of  micro  channel  and  Reynolds  number  on  methanol 
conversion  and  H2  production  rate  at  Tw  =  260  °C. 
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catalyst  loading  and  a  slower  velocity  than  micro  channel  (II). 
This  is  because  that  micro  channel  (II)  has  a  more  uniform  tem¬ 
perature  distribution  which  increases  the  chemical  reaction  rate. 
Comparison  of  Fig.  11  (a  and  b)  shows  that  the  CH3OH  consump¬ 
tion  increases  with  increasing  wall  temperature,  which  in  turn, 
increases  the  reaction  rate. 

5.  Conclusions 

A  three-dimensional  micro  channel  model  was  developed  to 
investigate  the  geometric  sizes  (aspect  ratios  and  micro  channel 
size)  and  thermo-fluid  parameters  (Reynolds  number  and  wall  tem¬ 
perature)  on  methanol  conversion  and  local  transport  phenomena 
in  the  micro  channel  of  a  plate  steam  methanol  micro-reformer.  The 
conclusions  are  as  follows: 


Fig.  10.  Effects  of  geometric  size  of  micro  channel  on  local  temperature  along  center 
line  of  micro  channel  at  Tw  =  200  °C  and  Tw  =  260  °C. 


it  can  be  shown  that  the  local  temperature  distribution  increases 
with  increasing  wall  temperature. 

The  effects  of  the  Reynolds  number  Re  on  the  methanol 
mole  fraction  distributions  of  two  micro  channel  sizes 
(I)  33  mm(L)  x  1.5  mm(  W)  x  0.75  mm(H)  and  (11)33  mm(L)  x 

0.429(W)mm  x  0.214  mm(H)  for  wall  temperatures  of  200  and 
260  °C  are  examined  in  Fig.  11.  The  results  show  that  a  larger 
CH3OH  mole  fraction  is  noted  for  a  case  with  a  higher  Reynolds 
number  Re.  Additionally,  it  is  clearly  observed  that  for  fixed  Re, 
micro  channel  (I)  shows  similar  methanol  distributions  to  those 
of  micro  channel  (II).  Therefore,  micro  channels  (I)  and  (II)  have 
similar  methanol  conversion,  but  micro  channel  (I)  requires  more 


Fig.  11.  Effects  of  geometric  size  of  micro  channel  and  Reynolds  number  on  local 
CH3OH  mole  fraction  along  center  line  of  micro  channel  at  (a)  Tw  =  200 °C  and  (b) 
rw  =  260  °C. 


1.  A  reduction  in  the  aspect  ratio  (y)  improves  FI2  production  rate 
and  methanol  conversion. 

2.  A  decrease  in  the  aspect  ratio  can  increase  methanol  consump¬ 
tion  and  the  hydrogen  mole  fraction  along  the  micro  channels. 
A  reduction  in  the  aspect  ratio  also  yields  a  much  more  uni¬ 
form  temperature  distribution,  while  also  increasing  the  steam 
reforming  reaction  rate  and  methanol  conversion. 

3.  The  better  methanol  conversion  and  the  H2  production  rate  are 
noted  for  the  smaller  aspect  ratio. 

4.  Smaller  micro  channels  have  much  more  uniform  temperature 
distributions,  which  in  turn,  result  in  a  decrease  in  the  ther¬ 
mal  entrance  lengths.  Therefore,  fuel  utilization  efficiency  is 
improved  for  a  small  micro  channel. 

5.  For  the  same  hydraulic  diameter,  a  reduced  Re  increases  the  reac¬ 
tant  gas  residence  time,  which  in  turn  increases  the  reaction  time 
and  significantly  increases  methanol  conversion. 
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